Fermion Superfluids of Non-Zero Orbital Angular Momentum near Resonance 
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We study the pairing of Fermi gases near the scattering resonance of the I 7^ partial wave. Using 
a model potential which reproduces the actual two-body low energy scattering amplitude, we have 
obtained an analytic solution of the gap equation. We show that the ground state of £ = 1 and 
£ = 3 superfluid are orbital ferromagnets with pairing wavefunctions Yn and Y32 respectively. For 
£ = 2, there is a degeneracy between Y22 and a "cyclic state". Dipole energy will orient the angular 
momentum axis. The gap function can be determined by the angular dependence of the momentum 
distribution of the fermions. 



The discovery of fermion superfluids near Feshbach res- 
onances is a major development in cold atoms physics 
Not only the long sought goal of achieving a fermion su- 
perfluid has been realized, the resulting condensate ex- 
hibits many interesting "universal" properties due to the 
interplay between unitarity scattering and Fermi statis- 
tics |2j. Even though the studies of these superfluids 
have just begun, another new and exciting direction has 
already been spun off. This is the physics of resonances 
with non-zero orbital angular momentum. 

While most experiments on Feshbach resonance fo- 
cus on s-wave scattering, there are many Feshbach reso- 
nances with non-zero momentum accessible by the same 
energy tuning method using external magnetic fields. Re- 
cently, Salomon's group at ENS has reported evidence of 
reversible production of p-wave molecules by sweeping a 
Fermi gas of 6 Li through a p- resonance 3] . It is therefore 
conceivable that p-wave or even higher angular momen- 
tum fermion superfluids can be realized in the future. 
Typically, the width of the resonance decreases with in- 
creasing angular momentum. One therefore expects that 
I 7^ superfluids will be harder to observe than their s- 
wave counterpart. Whether this can be achieved within 
current technology remains to be seen. The fact that the 
width of p-resonances can be resolved in recent experi- 
ments |3] is very encouraging. In any case, the difficulty 
is not an intrinsic one. We hope that the novel properties 
of the £ 7^ superfluid pointed out below will motivate 
searches for these new superfluids. 

In this paper, we consider I 7^ pairing near a scat- 
tering resonance where interaction between particles is 
strongest. Our goal is to determine the ground state 
structure, their signature, and the possible existence of 
universal behavior in these systems. As a first step, we 
shall focus at T = 0. In the case of s-wave resonance, 
it has been shown ,4s] that mean field theory is valid at 
T = 0, despite large fluctuation effects near T c . Using 
a model potential that reproduces the exact two-body 
scattering amplitude in vacuum, we have found analytic 
solutions of the BCS problem for all £ 7^ pairing at 
T = 0. Our findings are: (A) For p and /-wave pairing, 
the pairing states are Yn and Y32 respectively, which 
are orbital ferromagnets that break time reversal sym- 



metry and carry macroscopic angular momenta. For d- 
wave pairing, there is a degeneracy between Y22 and a 
so-called "cyclic" state. (B) The criterion for determin- 
ing the ground state structure of a £ 7^ superfluid is 
that the energy gap must have minimum angular fluctu- 
ation. (C) Unlike s-wave superfluids, where energy per 
particle has a universal form £ F (1 + (3), the properties of 
£ 7^ superfluids are not universal and are determined 
by the effective range r, of two body scattering. The 
energy per particle is proportional to £p\kFfi\ << £p , 
where k F is the Fermi wavevector. (D) Although dipo- 
lar energy is insignificant for superfuid pairing, it breaks 
the rotational symmetry and orients the angular momen- 
tum of the pair. (E) Experimentally, the nature of the 
pairing state can be easily revealed by the angular de- 
pendence of the momentum distribution of the fermions, 
whose orientation can be controlled by an external mag- 
netic field through the dipole interaction. These results 
are established below. 

To begin, we first point out a major difference between 
the £ 7^ pairing of atomic gases and the p-wave pairing 
of superfluid 3 He. In the latter case, the pairing interac- 
tion is rotationally invariant in spin space. The p-wave 
interactions between 3 He pairs (TT),(T4 + IT), (44) are 
identical. In the recent ENS experiment, three different 
p-resonances are found for the 6 Li pairs in spin states 
||^)> and \ = ^ = t)- When the pair ||^) is at 
resonance, the pairs \\\), and are not and their 

interactions can be ignored. In other words, particle in- 
teraction is highly anisotropic in the "pseudo-spin" space 
for atomic Fermi gases near resonance. 



Setting up the pairing problem: We consider a 
two component Fermi gas (denoted as \ and 4) with a 
Hamiltonian H — H D + V, H — Y^i 



H 2 k 2 /2M, M is the mass of the fermion, V 
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where Vq = J dre~ zcl ' T V(r) is the Fourier transform of 
the potential V(r) between unlike fermions, and f2 is the 
volume. Interactions between like fermions will be set 
to zero. If b is the range of the potential V(r), then for 
wavevector k such that kb << 1, the two body scattering 
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amplitude in the £-th partial wave channel is 

(kb) 2e 



/<(*) 



-aj 1 +r e k 2 /2 - i(fcfe) 2£ fc' 



(1) 



where a t and r t are the scattering length and effective 
range respectively. At resonance, ag diverges, while rg 
remains of microscopic size. In the case of a square well 
(V(r) = —\V\ for r < b, and otherwise), it is straightfor- 
ward to show that near resonance r e = — (21— l)!! 2 . 
Recall also that the bound state energy E b will appear 
as a pole in f(k) when k is continued analytically to the 
pure imaginary axis, k — > in. Near resonance, we have 



E,, 
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(2) 



It is clear that a bound state exists only when a £ r e < 0. 

The scattering amplitude fg(k) is related to V(r) 
through the T-matrix Tg(k, k'; E), 

M 

h{k) = -—^T t {k,k;2e k + iQ + ), (3) 
Awn 

where Tg(k,k';E) satisfies the integral equation 

dp p 2 V l (k,p)T i (p,k';E + ) 



T i (k,k';E+) = V e (k,k') + 



2tt 2 E + -h 2 p 2 /M 



(4) 

with E + = E + i0 + , where we have used the expansion 

Vfc-k = 4^£, V t {k,k') YL=-t Ytmfa Yf m (&). 

In standard BCS theory, the ground states is \G) = 

]lk(uw+Vual^aL k>l )\0), |uk| 2 + |«k| 2 = 1- The coherence 
factor u k and v k are determined by minimizing 

{H-fiN) =2^(e k -/i)K| 2 -0- 1 X)^-kU k , Uk ,< Uk) 

k k,k' 

(5) 

which gives » k = y/(l+ £jE k )/2, v k /u k = A k /(£ k + 

£k)> = V^k + l A k| 2 > C k = e k - l*>, and the energy gap 
A k satisfies 



A k = - ^ V k-pUk>V p = 



p 



Vk-pAp 
2E„ 



(6) 



The chemical potential fi is determined by the number 
constraint 



N = 1 ^ / e k - n 



(7) 



Since we are interested in the resonance physics of 
the £-th partial wave, we replace Vk'-k simply by its 
£-th angular momentum component, i.e. Vic'-k — > 
47r^(fc, fc') EL-/ Y em (k)Y; m (k>). ^ 

Even with this (single harmonic) simplification for 
Vic-ic', there is a serious complication which does not 
occur in s-wave pairing that eq.@ does not have a 



single harmonic solution, since nonlinearity will force A k 
to have all spherical harmonics. While one expects on 
physical grounds that only a few harmonics around £ are 
dominant, there is no simple way to extract such domi- 
nant piece and to calculate the dominant energy contri- 
bution. 

To eliminate this technical complication, we adopt the 
viewpoint that all microscopic potentials that produce 
the same low energy scattering amplitude will describe 
the same low energy physics for the system. We can 
therefore replace the actual Vk, k ' by a model potential 
that produces the same scattering amplitude f e , but al- 
lows a much easier solution for the gap equation. The 
most convenient model is a generalization of the separa- 
ble potential used by Nozieres and S. Schmitt-Rink |g, 



V e (k,k') = Xtw t {k)w e (k'), 



wg{k) = 



{k/k 



[l + ik/ko] 2 ]^ 1 )/ 21 



(8) 



(9) 



where k is a momentum cutoff. With eq.JSJ), eq.Q has 
the solution 

T e (k, k'; E+) = n{E + )wi{k)w t {k'), (10) 



1 1 1 ^ w e (k) 2 



te(E+) ^ k — 2e k 



(11) 



For £ 0, eq. jTI} has a low energy expansion Q, 
1 

k 

Y J Mk?5{E~2e k ). (12) 



_ 1 l^ wgjk) 2 E^ wgjk) 2 2 

tg{E + ) \e O 4- 2e k + Q ^ 4e 2 + U ^ " J 
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Substituting ea.ljTU)! and lfT2^) into eq.©, and not- 
ing that the last term in ea. l|12|l integrates to 
ikw^kf^T^/M)- 1 , we achieve the form eq.Q pro- 
vided A and k Q are related to the physical parameters ag 
and r„ as 
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More explicitly, we have 
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g 2)(£+l)n' 



(16) 



Near resonance, a t — > oo, the second term in r g in ea. 1(14(1 
and 1(15(1 can be ignored. 

Solution of gap equation: With the potential given 
by JHJ, eq.® has a solution 



Yi m (k), 



(17) 



where the coefficients {C m } satisfy the equations 
-C m /Xt = Anil- 1 £ p , m ,^(p) 2 YL(f>) YirrAP) c m ,/ 
(2E p ). Using ea. fH^ to express \t in terms of physical 
parameters a £ and r e , we recast this equation as 



Cm 

9e 



-in 



p.m' 



2E p 2e p 



The energy of the system, from eq.©, is 



|A k | 2 



2E 



fj,n. 



(18) 
(19) 



It is useful to separate the angular structure and mag- 
nitude of A k by writing 

A k = w e (k)C A(k), A(k) = ^(k), (20) 



where C 2 = £ m |C™| 2 ana a m = C m /C. Eg. (fT5|) can 
then be written as 



C 2 



/ 1 



Our goal is to solve from eq. ((21(1 and J7J) the quantities (C 
and A(k) and /i) as a function of (n, g e , and r^), and then 
determine which solution (C, A(k)) minimizes ea. (|19H . 

In the following, we shall present an analytic solution 
for ea. 1(21(1 and Q) near resonance. We have also solved 
these equations numerically for all gt. For all the angular 
momenta we have studied, the two results are identical 
within the region of gt where the analytic solution is valid. 
To derive the analytic solution, we write ea. (|21Jl and 
Q in dimcnsionlcss form by expressing all energies and 
wave- vectors in units of e ko and k Q , i.e. defining \i = ~p.t kol 
C = Ce ko , k = qk D , e k = q 2 e ko , A k = e ko w e (q)C A(q). 
Next, we assume that the solution of ea. (|21J) and 10 for 

2 

1^0 satisfy ~p, C < < 1 . We will verify later that this 
is indeed the case when k F r e \rg/a e \ « 1, With this as- 

sumption, we can expand eq. (|21|l and Q) in ~p and C . To 
the lowest order in these quantity, eq.JSlJ, J7J) and i(TT)|) 

(when expressed in terms of un-scaled variables) become 



4:irg g fl 



-y 



|A P I 2 

K 



i 



\ - 



|A P | 4 



(22) 



-V 



|A P | 2 
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(23) 



(24) 



In deriving eq. 124|) , we need to include second order terms 

„ 2 4 

~p , ~pC , C due to cancelation of lower order terms. 
Note that these expansions will not work for s-wave since 
the sums in ea. (|22Jl to 1)23(1 are infrared divergent. 
Combining ea. 1(22(1 and 1(23(1. we have 



n\i = n/j, - - , nfi _ 

47rg £ J I 



o 



|A r 



4e3 



(25) 



Using ea. ((23(l and ea. ((14(l near resonance, (and noting 
that the l/(k ae) term in ea. l(14|) can be ignored near 
resonance), we find two equivalent forms for C, 



C 2 



nE b 
2 : 



C 



S(k Q b) e 



4arg e a e r e M 2 ' ^3 a/M 

Comparing ea. l(25|) and 126(1 . we have 
H = f i + EJ2. 



(26) 



(27) 



The explicit form of \i a can be obtained by evaluating 
the integral in ea. ((25(l and using the expression of r f in 
eg. ((15(1 near resonance, which is 



Mo = 7 



/ 2^ 2 I 2 



V3(fc 6) 2£ 



d P |A(p)r 



(28) 

where I„'s are given in ea. lfTH|) . Ea. (|2l))l to (f^ to- 
gether with eq.® give C and fi as a function of n, a £ 
and r e . From these equations, it is easy to show that 

2 

~p,C ~ I^f^I 3 - Since |fc F r£| << 1, our initial assump- 

2 

tion JL,G << 1 is valid. Note that the structure of the 
gap 7 shows up only in fi and not in C . 

Finally, using eq.J31|), HU, anQ (EH > we obtain the 
energy density as a function of n, and r^, 



(H)/Q = nu>- nn /2 = (n/2)(^ + E b ). 



(29) 



Eq.(|25I), pgjl. and (01 imply that (i) the ground state 
has minimum angular fluctuation in the gap, i.e. 7. (ii) 
Unlike the s-wave case where \i is of order £ F near res- 
onance, ea. 1(28(1 shows that the chemical potential /i for 
t =^ at resonance is greatly reduced from £p, by a non- 
universal factor |/c F r^|, reflecting a stronger interaction 
energy then the s-wave case. This is due to the fact that 
the energy of the bound state E b for t ^ grows much 
faster than that of s-wave away from resonance, since 
E b oc l/(a e r e ) for I 7^ whereas E b cx -l/a§ for £ — 0. 
The effect of this larger interaction energy also shows up 
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FIG. 1: Chemical potential (fi) for I = 1 resonance as a 
function of — 1/fcFdi. We have taken n = —36, k given by 
eg. 1151 . and k F r t = —0.09. The linear portion agrees exactly 
with eg. 1271 and (1281) . The switching to /i = £ F at the atomic 
side is very sudden and is controlled by the smallness of k F r 1 . 
It takes place around Eb/2 ~ £ F , or (kpd^ 1 ~ k F r 1 . 

at high temperatures. However, the thermodynamics in 
that regime is universal 

Ground state structure: To find the mini- 
mum of 7, we shall use the rectangular representa- 
tion of spherical harmonics to write ^ m a m Y« m (k) = 
^-ii,i2,---,itK 1 K 2 ---K e ^ where A is symmetric in all 
its indices, and vanishes upon contraction of any two 
pairs of indices. For example, the sum ^2 m a m Yi m (^) 
is A • k, A a ^ b k a k b , A a ^ c k a k b k Cl for £ = 1,2, and 3 re- 
spectively, with TiA = for 1=2 and J2 a Aa,a,c = 
for I = 3. We then have j e=1 cx 2(A* • A) 2 + |A • A| 2 , 
7 f=2 oc 2(TrA t A) 2 + |TrA 2 | 2 , and a similar but lengthier 
expression for £ = 3. 

For £ = 1, the minimum occurs at A 2 = 0. This means 
A cx x + iy or any rotation of it, which is A(k) = ^ii(k) 
along an arbitrary angular momentum quantization axis. 
The superfluid is an "orbital ferromagnet" since all pairs 
are in the same Yn state. This is remarkable for it im- 
plies that by sweeping across the Feshbach resonance, 
a superfluid with macroscopic angular momentum and 
broken time reversal symmetry will result. 

The problems of minimizing 7 for t = 2,3 were 
solved by N.D. Mermin in the context of superfluid 
3 He 0, El- For £ = 2, there is an accidental de- 
generacy. Both A(k) = l22(k) and the "cyclic" state 
A(k) cx fc 2 + e 27 ™/ 3 fc 2 + e 4 ™/ 3 *; 2 minimize lt=2 0. This 
degeneracy can be resolved in higher order in | A; t~2 | and 
will be discussed elsewhere. For £ = 3, the ground state 
is A(k) = Y32(k) along an arbitrary direction [Til . This 
state is also an orbital ferromagnet, even though it is not 
of maximum angular momentum state. At present, there 
are no exact solutions for £ > 4. 

Numerical Results: Note that although fj, ~ 



S F \k F r e \ near resonance (eq.((35Jl), it has to recover to 
(i ~ £ F on the atomic side of the resonance whether a t 
is negative and small. We have solved eq.J2U and (Q) 
numerically and have shown this is the case. (See figure 
1). Our numerical results are in exact agreement with 
ca. (|2*7Jl and (|2*S|) near resonance. 

The effect of dipolar energy: Dipolar energy V D 
breaks rotational symmetry in real space. Since electron 
spins are polarized by the external magnetic field B, we 
have V D = i/£/(r-r')Vi(r)^(r')^(r')V'a(r), U(v) = 
H 2 B {1 — 3(B • f) 2 )/r 3 , where \ib is the electron Bohr mag- 
neton. Since dipolar energy per particle is fi|n, and since 
H%n/n = {e 2 k /m e c 2 ) (M/m e ) (3tt 2 /2) - lfr 4 for k a ~ 
10 cm -1 , dipolar energy is not strong enough to affect 
the gap structure. On the other hand, it can orient the 
pairing state. A straightforward calculation shows that 
(Vb) = (2^ 2 3 /f7)E k:k ,(l/3-(B.q) 2 )A*,A k /(4 J B kJ B k ,), 
where q = k' — k. In the case £ = 1, A k cx A • k, we 
have (Vr>) oc — |B • A| 2 . Since A cx x + iy, we have 
(Vd) cx +(B • z) 2 . The angular momentum of the pair 
will lie in the plane perpendicular to B. 

Signature of the I 7^ superfluid: Eci. (|23|l shows 
that the momentum distribution n p of the fermions is 
|A p | 2 /(2e 2 ) cx |A(p)| 2 . A measurement of the angular 
dependence of n p therefore gives |A(p)| 2 directly. 

We have thus established results (A) to (E) mentioned 
in the Introduction. This work is supported by NASA 
GRANT- NAG8- 1 765 and NSF Grant DMR-0109255. 
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